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Introduction 



Abstract 

Brans-Dicke scalar-tensor theory provides a conformally coupling of the scalar field with gravity in Einstein's frame. 
" " j " This model is equivalent to an interacting quintessence in which dark matter is coupled to dark energy. This provides a 
I . natural mechanism to alleviate the coincidence problem. We investigate the dynamics of this model and show that it leads 
' to comparable dark energy and dark matter densities today. 

cm : 

&; 

^According to Einstein's theory of relativity, the observed late time acceleration [JJ of the universe can be successfully explained 
Qy introducing the dark energy component [5] on the right hand side of the Einstein's equations. There are two main scenarios 
used to explain the nature of dark energy: a cosmological constant A[3j or a scalar field, usually called quintessence [J. The 
i— ec>smological model that consists of a mixture of cosmological constant and cold dark matter is called ACDM and the other 
^ po ssibility leads to QCDM cosmological model. Although ACDM model opens the possibility to fit it with the observational 
|Uata, two questions arise: first, why the value of the cosmological constant is so small, and second, why it's energy density 
^gins to dominate now. This is known as coincidence problem [5]. 
Quintessential models of dark energy also provide an adequate fit to the cosmological data. Moreover the existance of 
stable tracking solutions of the quintessence field shows the importance of this model in the investigation of the fine tuning 
' ^problem. But it can not explain the present ratio of energy densities in the universe, i.e. the coincidence problem [6]. One 
Q^iossible way to make a model to overcome the coincidence and fine tuning problems, is to introduce a proper coupling term 
O^ietween the dark energy and dark matter. The form of the interaction term is motivated from the phenomenological consid- 
eration or by requiring that the ratio of dark energy to dark matter densities, r = p v /p Cl be stable under perturbations [BJ [7]- 
CM ' Introducing the interaction in QCDM model, can lead to a finite r at present or causes this ratio varies slowly now, 
QYhich is called soft coincidence [8]. Moreover it is required to have a smooth transition from matter dominated universe to 
^^h'e accelerated expanding universe. Denoting the coupling term by Q in the conservation equations of the quintessence and 
^T*osmic dark fluid, one has: 

Pc + 3Hp c = Q (1) 
• S ! p v + 3H(l + u)p v = -Q (2) 

which p c , p v = \<p 2 + V((p) and p v = \Cp 2 — V(tp) arc the energy density of cosmic dark fluid, the energy density and 
pressure density of the scalar field respectively, lo = p v j p v is the equation of state parameter for the scalar field and a dot 
over any quantity denotes derivation with respect to the cosmic time. 

A variety of interaction terms have been suggested in the literature [S [7l [8l [9l [TU l [TT ] [T2 l [13 l [T4 1 [T5 l rT6 | rTT l rTH j ri9 l [20]. 
One of the most important ones is Q = \\Hp c [5J [TUJ HU Ell- Inserting it into the above equations, together with the 
Friedmann equations, yields a monotonic decreasing function of the scale factor for r which changes slowly at the present 
time, i.e. r/r\ t <C Hq [lOj . Furthermore assuming an accelerated universe and negative r, together with satisfying the soft 
coincidence criterion, the functions H{a) 1 V(tp) and the allowed range of Ai are determined. More general interaction term 
Q = \\Hpc + XiHp^p is investigated in [9j [10]. This yields to the evolution from a stable value r + at early times to an 
unstable value r_ at late times [TUJ. Also the constants Ai and A2 are constrained by the evolution of the background energy 
density. 
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It must be noted that the observational data are unable to discriminate between this interaction and the previous one. 
Also these two models are investigated in a more general manner in which the transfer rate of energy densities is determined 
by the properties of dark sector and in general is different from the expansion rate of the universe[13]. 

In [12], it is assumed that the scaling solutions (having a constant r at all times) exist for the generic form of a non- 
minimally scalar-field Lagrangian. This leads to the coupling Q = f3p c tp in the Einstein's frame. It can be shown |15| that 
the proper sequence of cosmological epochs in the form of radiation-dominated, (^-matter-dominated and accelerated epoch 
is followed in this model. Also the cosmological consequences of this kind of coupling is studied in P31 [151 EH HZ] assuming 
an exponential potential for the quintessence field. 

The simplest model which can provide this particular interaction is a non minimally coupled gravity in which the scalar 
field couples minimally to gravity but couples conformally to the matter fields [15l [18]. This can be obtained from the 
following action |21) : 

Se — J d A x^[R - i(V(^) 2 - V(cp) + C c (e- 2 ^g^ v ) + C b ( 9fiv )}. (3) 

It is obvious that this model is conformally related to a non minimal gravity in the form of Brans-Dicke type lagrangian. 
The cosmological properties of this model is investigated in the literature. In [22] an observational analysis of non-minimally 
coupled theory and it's differences from Genaral relativity is presented. Also the authors of [23] discuss the attractive aspect 
of non minimal quintessence and show that after tracking era, the quintessence dominated epoch is started. The particle 
production of this model in inflationary epoch is studied by Sahni et. al. |24j . Their analysis leads to a spatially flat universe 
consisting of dark energy (in the form of quanta of the non-minimally coupled scalar field created during inflation) and dark 
matter. In this way it presents a new way of generating dark energy quantum mechanically [35]. 

For the latter coupling, more general form of the potential and allowing /3 to be a function (3(ip) can be found in [19] . 

In summury, the coupling of dark energy to dark matter shows that the dark energy can affect the expansion history 
of the universe and the change of the growth rate of matter perturbations in the structure formation stage. Therefore this 
makes tight bounds on the strength of coupling terms [3D] . Also in [TT] , the different forms of coupling and their aspects have 
been investigated and compared with the observational data, the result indicates that the parameter space, which includes 
the arbitrary coupling constants and the equation of state parameter, is greatly limited. 

This work rests on the latter coupling pointed above. Moreover we use new variables, proposed in |10j . to simplify the 
analysis of the critical points of the non-minimal scalar field model. 

The outline of this paper is as follows. In the next section we shall study the critical points of the autonomous system 
of equations associated with the non-minimally coupled quintessence model. Then in the third section we disscus the 
dependence of dark energy-dark matter ratio on the redshift. Finally in the last section we shall give some conclusions. 



2 The model and its critical points 

Here we consider the quintessence field decays into the pressureless CDM according to the latter kind of coupling term 
discussed in the introduction, i.e. Q = Pp c <p. We also include baryons' density as an independent conserved component such 
that: 

p b + 3H Pb = (4) 

The conservation equations (TT]), ^ and Q can be obtained from the action It has to be noted that similar models are 
disscussed in |15| (in which baryons plus CDM couples to dark energy component) and in |16j (neglecting baryons). Both 
these works assume an exponential potential for dark energy Here we shall use a more direct way to investigate the model by 
assuming that the equation of state parameter, uj, is constant and then absorb it in the governing equations via introducing 
new variables, and we have not consider any assumption about the special form of the potential. 

Now, consider a spatially flat FRW universe occupied by the above three components. The cosmological equations of 
motion are (we have chosen the unit in which 8ttG = c = 1): 

H 2 = \{p v + Pb + Pc) (5) 
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Table 1: Critical points of the model and their properties. 
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Defining the dimcnsionlcss variables: 



H = --W(l + oj)+p b + p c ] (6) 



/V Pc Q 

• T= 3iP y= 3iP U= 3i^ (7) 
and following the same method as [TU], the conservation relations of dark sector, equations ((T|) and © can be written as: 

x = 3x(l — x) — u, (8) 
y' = — 3xy + u (9) 

in which u = —u/lo and a prime denotes derivation with respect to A = —win a (we assume the current value of the scale 
factor is one, ap = 1). Using the relation ip 2 = (1 + iS)p v , one obtains: 



u = ±Pyy/3(l+Lj)x (10) 

where f3 = /3/\lj\ and the positive (negative) sign indicates that ip and oj have the opposite (same) signs. Here we consider /3 to 
be a positive constant. Inserting (jTUJ) in ((5J and ((9]), the critical points of the above dynamical system together their stability 
along the dynamical state of the universe are analyzed by means of eigenvalues method and the results are summarized in 
table (UJ. 

The critical points are labeled by A, B and C and the type / (II) solution corresponds to the opposite (same) sign of <p 
and u>. The first critical point, A, is (x = 0, y = yo) which results in u = and thus there is neither scalar field nor interaction. 
This point represents the baryon-cold dark matter domination at early times, with no contribution from dark energy. For 
this point, the secular equation of the autonomous system (|5J and (J5J is singular. A detailed study of the equations near 
this critical point shows that the system diverges exponentialy from this critical point, and thus it is an unstable point. 

For the second critical point, B, we get (x = 1, y = 0) and again there isn't any interaction, u = 0. This corresponds to 
the pure dark energy dominated era which can be accelerating if uj < —1/3. Therefore depending on the initial conditions, 
approaching this point, the final state of the universe can be a pure dark energy dominated one in this model. 

Finally the third point, C, corresponds to the dark energy-cold dark matter dominated universe. It is only of type I, so 
that the sign of uj and ip is different. It can represent the late time universe as the previous one. For the latter points, B 
and C, one of the eigenvalues is always negative. The second eigenvalue for model B Type II is negative, while for models B 
Type I and C an appropriate range for (3 can be chosen so that the eigenvalue becomes negative. Therefore model B Type 
II is stable, while models B Type I and C are stable or saddle depending upon the value of the coupling constant, f3. On the 
other hand the critical points should be such that < x < 1 and < y < 1 which for the point C is satisfied provided that 
(3 < J j2—\ui\. Therefore the point C is a saddle point. We sec that there are only two critical points that admit accelerated 

solutions, B and C. These are in agreement with the solutions a and bu of Amcndola in [15] . The point C is the only 
accelerated solution in which the ratio of matter and dark energy densities is constant and therefore is a scaling solution. The 
existence of this critical point can provide the same order of energy densities for both interacting components independently 
of the initial conditions. This point is also in agreement with the solution D of Bohemer et.al in [TB] (See table I of |16|). 
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Restricting ourselves to the positive value of Q, which means tp > 0, corresponds to the transfer of energy from quintessence 
to the dark matter, and this is compatible with the second law of thermodynamics [26J. In this case, allowing the observed 
acceleration expansion of the universe, we have to select the type / solution from any critical point. 



3 Dark energy— Dark matter ratio 

In what follows, we shall discuss the evolution of the ratio of dark energy to dark matter densities, r. To do this, from 
equations (JH|) and we obtain [ID] : 

r' (r + 1) 

- = 3- K > u, (11 

r ry 

With the interaction term given by (1101) . the above equation becomes: 

- = 3toV^(1 + -), (12) 
r r 

where a = /3^/3(l + lj). To solve this equation, one needs the present value of r as the initial condition. However as the 
authors of point out it depends on the present time value of Hubble parameter and on the other hand there aren't any 
model independent observational data for r(H) and r(z). Thus these authors propose to determine the dependence of the 
Hubble parameter in terms of redshift, z, to compare with observational data which arc model independent |27j . To do this, 
using relation ©, equation ([T2"j) becomes: 

dr 2o;(3rT5Vi(r + l)) 

dH 3H(l+ux) ' 1 j 

Now combining equation (|12[) with the expression H = —dz/[(l + z)dt], yields: 

dr uj(3r =p a.y/x(r + 1)) 
dz ~ 0- + z) 

Dividing equation (TT4"]) by (jTU)) . gives: 

dH _ 3H(l+ux) 

dz ~ 2(1 + z) ' ( ' 

This relation shows that H(z) depends on the interaction term, Q, through the dark energy density in this model, see relation 
©• 

Although we have not an exact solution of equations (ITil) and (|15|) in general, the behaviour of r can be studied for 
extreme limits of the early and the late times. We shall show that provided that uj and ft are chosen appropriately, we have 
solutions that predict the correct behaviour of dark-energy-dark-matter ratio. 

At early times (matter dominated era), for which x — > 0, equation (fT5)) solves to H ~ (1 + z) 3 / 2 , and equation (fT4|) reads 

as: 

dr 3a;?' 

~dz = T+~z ^ 

with the solution r ~ (1 + z) 3uj . 

On the other hand, at late time, the baryonic matter density is ignorable and thus x + y — > 1, equation (|14[) reads as: 

dr 3ur auj 



(14) 



dz 1 + z 1 + z 

which has the following exact solution: 



Tr—y/r(l + r) (17) 



l + z=( yffT 1 ^ (18) 
1 ^i(r )F 2 (r )y V ^ 
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where: 

Fi(r) = {l + 2r + 2^/r(l + r)) ±a (19) 

and 

^a(r) = (a (-6^(1 + r) + a) + r (9 + a 2 )) (20) 

and ro is the current value of r. 

Now we shall set some assumptions on the asymptotic value and slope of dark energy-dark matter ratio. The late time 
value = R should be such that: 

R > ro (21) 

where ro is the present value of r. Moreover to have a nearly constant dark energy-dark matter ratio at late times we have 
the following restriction: 

£ « 1 (22) 

aZ r=R,z=Z 

where Z tends to —1. 

Using equation (fl7|) . the conditions (|2~T|) and (|22|) lead to: 



and 

w « l + , Z (24) 

3^-5^(1 + 7?) 

The evolution of dark energy-dark matter ratio with respect to the redshift is shown in figures ([1} and © . In figure ([T]) 
graphs for different values of /? and a fixed uj are drawn while in figure ^ graphs for different values of lo and a fixed /3 are 
considered. We see that this model leads to a constant unstable r at late times and to a lower constant stable r for early 
times. This shows a similar behaviour to what is obtained in [TUl H7J for two broad classes of interacting models of dark 
energy. 

In figures ([3]), (|4]) we have plotted the phase space diagrams of type / model (equations (JSJ), © and (|10[) V Figure ([5]) is 
the corresponding plot for type II model. 

In figure ([3]) the stability condition of BI point is satisfied and thus the point CI is a saddle point as it is expected from 
table 1. Changing the value of f3 such that the stability condition of BI point is not satisfied leads to figure Q and one can 
see that the critical point CI is out of acceptable range of x and y. 

Finally as it can be seen in figure ([5]), for a type II solution the point BII is an attractor point. 

A question that arises here is that which kind of potential leads to the derived solution. In order to reconstruct the 
potential, one can use the relations V = (1 — uj)p ip /2, (p = + uffhp and equation ^ to obtain: 

dV 1-uj 



dip 2 

p v in the above equation can be substituted by integration equation ([IJ leading to|14j: 

3 

«0 



'm^{l + Lo)p^ + Pp c \ (25) 



Pc = Pco ( — ) cxp (/% - (p )) (26) 



where a subscript denotes the current values of quantities. 

At late times, z ~ — 1, p v / 'p c ~ R and using the equations ([T]) and © we get: 
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Figure 1: Dark energy-Dark matter ratio as a function of redshift. uj is chosen to be —0.5 and (5 equals to 1.14, 1.19 and 
1.16 for dashed, dotted and thick lines respectively, r at present time is chosen to be 7/2.5. 

Inserting the above equation in the equation^ leads to: 



For early times, p v is small and thus ip ~ and tp ~ po and therefore the above potential can be expanded up to first 
order. Using this approximate potential one can check that the Friedmann equations and the equations governing densities 
are satisfied. 

4 Conclusion 

In this paper we have investigated non-minimal quintessence model of dark energy. Following |10j , we have defined a new 
variable, N = — wlna, which reduces the parameter space into a two dimmensional space. Our results are in agreement with 
those of [ini US]. We found some range of the interaction parameter f3 given by the stability criteria of the critical points, 
and we presented an exact solution of dark-energy to dark-matter ratio for the early and late times. 




(28) 



and thus equation (j25j) can be integrated to: 




(29) 
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Figure 2: Dark energy-Dark matter ratio as a function of redshift. (3 is chosen to be 1.16 and ui equals to —0.49, —0.51 and 
—0.5 for dashed, dotted and thick lines respectively, r at present time is chosen to be 7/2.5. 

According to the early time behaviour (r ~ (1 + z) 3 ") the model behaves like ACDM model. On the other hand it can 
provide a finite dark-energy to dark-matter ratio at late times and thus it can alleviate the coincidence problem. In fact it 
is shown here that provided the strength of the interaction (/3) is chosen appropriate, this interacting model can predict the 
expected behaviour for the dark-energy to dark-matter ratio with a constant oj. 
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